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Abstract. We present a detailed evaluation of rj, the critical exponent 
corresponding to the electron anomalous dimension, at 0(1/Nj) in a large 
flavour expansion of QED in arbitrary dimensions in the Landau gauge. The 
method involves solving the skeleton Dyson equations with dressed propa- 
gators in the critical region of the theory. Various techniques to compute 
massless two loop Feynman diagrams, which are of independent interest, are 
also given. 
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1 Introduction. 



The important properties of a renormalizable quantum field theory are con- 
tained within the renormalization group equation from which one can, for 
instance, determine how Green's functions depend on the renormalization 
scale. Central to this equation are the (3- function and j(g), the anoma- 
lous dimension of the basic fields. Ordinarily one computes these within 
conventional perturbation theory, in a particular renormalization scheme, 
and hence determines the properties of the Green's functions to the same 
order. One difficulty with perturbation theory is that computations at suc- 
cessive higher orders become exceedingly tedious due in part to the increase 
in number of graphs to be considered and the complexity of the integrals 
which appear. Thus it is not possible to determine the higher order con- 
tributions to a renormalization group analysis with ease. One approach to 
alleviate this difficulty is to examine models in an approximation different 
from conventional perturbation theory, such as the large N expansion. In 
this approach the quantity gN, where g is the perturbative coupling con- 
stant, is held fixed as N — ► oo so that one remains in the perturbative 
regime of the model. Theories which admit such an expansion are those 
with an internal symmetry and included in this class is quantum electro- 
dynamics, (QED), with Nf flavours of electrons which we will examine in 
detail in this paper. From the graphical point of view the large N expansion 
is a reordering of perturbation theory such that chains of bubble graphs are 
summed first. As in conventional perturbation theory one can renormalize 
Green's functions and extract the pole structure from which one can deduce 
the large N approximation to /3(g) and 7(5). However, one will also run 
into the same difficulties as perturbation theory, such as the appearance of 
intractable integrals, which will occur at next to leading order. 

This problem is overcome by an alternative approach developed in ||, g] 
for the O(N) bosonic a model, which involves examining the theory at 
the d-dimensional critical point, defined as the non-trivial zero of the (3- 
function, where the theory is finite and massless. Moreover, the fields also 
obey asymptotic scaling, pL where the propagator will take the simple con- 
formal structure l/(x 2 ) a in coordinate space, with a its critical exponent. In 
particular, one examines the skeleton Dyson equations of the theory, which 
are valid at the critical point, and derives a critical point consistency equa- 
tion which can be solved within the large N expansion for the anomalous 
dimension, jl], g]. Since the anomalous dimension exponent, rj, of, say, a 
bosonic field is related to 7(5) via r\ = (d — 2) + j(g c ), where g c is the 
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critical coupling, then one can deduce the coefficients of 7(5) to all orders 
in perturbation theory within the particular large N approximation made. 
This and similar scaling relations emerge because the renormalization group 
equation takes a simplified form at criticality. Further, the absence of a 
mass for the basic fields means that one can quite straightforwardly probe 
the model beyond the leading order of the conventional large N approach. 

Recently this technique was applied to QED with Nf flavours of elec- 
trons in ||. There 77 was computed at leading order in the Landau gauge 
and so it is the purpose of this paper to present a detailed evaluation of the 
0(1/Nj) corrections to the result of j|] using some of the techniques intro- 
duced in earlier works, |IJ ||], f|, [|] , as well as developing others for the specific 
case in hand. Since the calculation of 77 is in arbitrary dimensions, we will 
therefore not only provide additional coefficients of j(g) in four dimensional 
perturbation theory in the MS scheme but also 0(1/ N 2 ) corrections for the 
three dimensional model, which is perturbatively super-renormalizable and 
currently of interest in various problems. As far as we are aware the only 
previous 0(1/ Nf) calculation in QED was carried out in ||, where 77 was 
computed in the Feynman gauge but precisely in three dimensions. Unlike 
the arbitrary dimension analytic result we give here, whose derivation is 
predominantly algebraic, the calculation of Q was partly carried out nu- 
merically. Moreover, the techniques developed here for four dimensional 
QED, will be very important when other more physically consistent theories 
like QCD are solved within this critical point large Nf formalism. Although 
QED is not a consistent theory in isolation due to the occurrence of the 
Landau pole at large values of the coupling this does not prevent us from 
accessing the perturbative region of QED in large Nf as then g is small. 

Earlier large Nf analysis of QED was carried out in [8| where the pole 
structure of the electron self energy and photon electron vertex were deter- 
mined by explicitly carrying out the large Nf bubble sum in the Landau 
gauge. The 0(1/ Nf) corrections to the /3-function were determined as well 
as the renormalization group function corresponding to the dependence of 
the renormalized mass with the renormalization scale both in MS, ji], H|, 
though the wave function renormalization was not studied. To go beyond 
this leading order by explicitly computing the next to leading order correc- 
tions to the bubble sum would be very involved and thus it is appropriate to 
follow the more efficient and elegant methods of [jl], ||] to compute 0(1/ Nf) 
corrections since it turns out that there are only two 2-loop corrections to 
consider. Finally, we refer the interested reader to previous 0(1 /N 2 ) calcu- 
lations in other models, ie pi 0, B, [|, M , since they will very much serve as 
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a basis for the calculation presented here. For instance, since QED involves 
fermions we will use several results from the much more straightforward 
0(1/N 2 ) calculation of r\ in the O(N) Gross Neveu model, Q, which is 
also fermionic. We note that the large TV" exponents which have been de- 
rived in these other models have all been shown to be in agreement with 
the appropriate renormalization group functions to as many orders as they 
have been calculated within explicit perturbation theory using dimensional 
regularization in the MS scheme. 

The paper is organised as follows. In section 2, we introduce the nec- 
essary notation and formalism we will use and review the leading order 
analysis of ||, but in coordinate space here. The formal corrections to the 
leading order consistency equation for rj are discussed in section 3, where 
we also deduce the renormalization required to analyse the theory in the 
critical region. Section 4 deals with some of the calculational techniques 
required and difficulties encountered, in explicitly computing the 0(1/Nj) 
correction graphs. We discuss their application to the photon self energy in 
section 5, as well as providing some detail on the calculation of the compo- 
nent bosonic graphs which arise and further techniques for their evaluation. 
A similar discussion for the corrections to the electron self energy are given 
in section 6, whilst we provide the main result of our efforts in section 7, 
where r\ is deduced at 0(1/N?). Our conclusions are contained in section 8. 
Finally, we provide various appendices which either list useful identities for 
general massless two loop Feynman diagrams or contain a library of the more 
involved basic two loop component graphs which occurred in our analysis. 



2 Preliminaries. 

We begin by introducing the formalism we will use in our analysis as well 
as recalling the important features of QED we require. To fix notation and 
conventions we will calculate with the (massless) QED lagrangian, 

L = W + A^W - %P - (2-1) 

where unlike the usual formulation we have rescaled the coupling constant, 
e, into the definition of the kinetic term of the C/(l) gauge field so that 
the dimensional analysis of each term in (2.1) is completely analogous to the 
earlier critical point treatment of other models, g]. Also, we have set 
= dp_A v — d y A^ and the field vp l , 1 <i < Nf, corresponds to the electron with 
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Nf flavours. The parameter, b, of (2.1) denotes the conventional parameter 
which appears in a covariant gauge fixing term. 

First, we recall the renormalization group functions, [3(g) and 7(5), of 
QED to the perturbative orders in which they are known in d = 4 — 2e 
dimensions e < small, where we set g = (e/2ir) 2 . Recently, /3(g) has been 
computed to four loops in MS in [01]] and is 



2N f „ N f o Nf(22N f + 9) 4 
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where we use the conventions of || |Tl|] and g is the dimensionless coupling 
constant in d dimensions. It is important to note that (2.2) is what is deter- 
mined as the /3-function in (4 — 2e)-dimensions in MS using dimensional reg- 
ularization prior to setting e = and that the coefficients are d-independent. 
(If they were d-dependent then one would not be using MS.) It is easy to 
see that when d < 4 (2.2) indicates the existence of a non-trivial zero, g c , 
of (3(g) which corresponds to a phase transition. To deduce the perturba- 
tive coefficients of the (4 — 2e)-dimensional functions of the renormalization 
group equation from the large Nf exponents we compute later, the location 
of the critical point is given from (2.2) by 




In (2.2), the coefficients of the large Nf term at each order in g agree with the 
explicit large Nf MS renormalization of the /3-function of QED carried out in 
[fL2f using dimensional regularization in d = 4 — 2e dimensions. Further, the 
critical exponent \jv = — f3'(g c ), which in effect carries information on the 
perturbative /3-function has recently been calculated using the formalism 
of |l|, ^] by examining (2.1) at the d-dimensional critical point, |13| , and 



is in exact agreement with the results of [12, [L0|]. The electron anomalous 
dimension, 7(g), is also known to several orders by restricting the QCD 
results of |l4|, to the abelian case and is 

b (4iV> + 3) 9 (40/V? + 54N f + 27) , 
1(9) = \g ~ K \l V + " V + 0(g±) (2.4) 



in a general covariant gauge. In fact it has been proved in [jl6| that the 
only 6-dependence which appears in the electron anomalous dimension is 
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at one loop and therefore the coefficients of the two loop and higher terms 
of (2.4) are gauge independent. Our remarks concerning the lack of d- 
dependence of the coefficients of the powers of g in (2.2) hold equally for this 
function since it was also computed in (4 — 2e)-dimensions using dimensional 
regularization with MS. In an earlier work, Q, we computed r\ = j(g c ) in 
the Landau gauge, 6 = 0, and deduced the 0(1/Nf) terms of 7(5) in the large 
Nf expansion to compare with (2.4). These were in total agreement with 
the coefficients appearing in (2.4) to 0(g 3 ) and subsequently we deduced all 
the higher order 0(1/Nf) coefficients. So, for example, in this gauge, ||, 

(35iV 3 \ 
—JL + hlN f + ClNf + d l \g 4 + 0(g 5 ) (2.5) 

where 61, c\ and d\ are unknown. By calculating 772, we will be able to de- 
termine the next to leading order coefficients of 7(5). Essentially, expanding 
7) at a particular order in 1/Nf in powers of e = 2 — ^d, the coefficients of 
e n in r\ are related to the nth order coefficient of 7(3) since g c ~ 3e/Nf from 
(2.3). Moreover, since the first non-trivial terms of the expansion of both rj 
at 0{1/Nf) and 0(1/Nj?) in powers of e begin at 0(e 2 ) they will produce in- 
formation to all orders on the gauge independent part of the renormalization 
group function. So in this sense we will regard 77 as gauge independent. 

Having recalled the perturbative structure of the field theory in d = 4 
— 2e dimensions and the equivalence of the leading order 1/Nf exponents 
already computed with the technique of 0, |2| for QED we now introduce 
the formalism of the method we use. First, we note the consequences for 
the Green's functions as a result of the existence of a non-trivial fixed point 
in d ^ 4 dimensions. From a statistical physics point of view, near a crit- 
ical point physical quantities obey simple power law behaviour where the 
power or critical exponent depends purely on the dimension of spacetime 
and the parameters corresponding to any internal symmetry by the univer- 
sality principle. (See, for example, fl7|| .) From the continuous field theory 
point of view, one is dealing with fields which are conformal at criticality 
and thus do not involve any mass. Therefore, we take the following coordi- 
nate space forms for the asymptotic scaling forms of the propagators of the 
fields of (2.1), which are consistent with Lorentz symmetry, near criticality 
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as 



I! §01 



ip{x) 




a 



B 



(2.6) 



{x 2 y 



(2/x - 2/3 - 1) x 2 



as x 2 — > 0. The point is that in this limit, or equivalently A; 2 — > oo, this struc- 
ture is the dominant part of the propagator which governs the renormaliza- 
tion and is therefore the part which is relevant for deducing the anomalous 
parts of the critical exponents a and (3 which are equivalent to the critical 
renormalization group functions. We have chosen to work in the Landau 
gauge for the following reason. Perturbatively, if one used any (covariant) 
gauge, other than the Landau gauge, then the gauge parameter gets renor- 
malized and thus the gauge changes. Therefore, since the large Nf expan- 
sion is a reordering of perturbation theory such that chains of bubbles are 
summed first, it is thus important that one works in a gauge, ie the Landau 
gauge, which is unaffected by renormalization effects, ||]. Whilst A^x) 
does not appear to take the usual form for the Landau gauge it is easy to 
transform (2.6) to momentum space using a Fourier transform to observe 
that its structure will be proportional to P(j, u (k) = (r]^ v — k^k u /k 2 ) where k 
is the momentum conjugate to x. For completeness and to fix conventions, 
we note that the Fourier transform we use is Q , 



where we have set a(a) = r(/z — a)/T(a) and the dimension of spacetime, 
d, to be d = 2fi. For the most part of this paper, we will work in coordinate 
space though it is straightforward to map from one space to the other via 
(2.7) since the fields are massless. The quantities A and B in (2.6) are the 
amplitudes of each field and are independent of x whilst a and (3 are the 
exponents of the respective fields. They are related to the exponents we will 
calculate, via, fH, 



where rj is the electron anomalous dimension and x 1S the anomalous di- 
mension of the electron photon vertex of (2.1). Both the latter quantities 
depend only on \x and Nf and are 0(1/Nf) within the large Nf expansion 
and will be calculated to 0(1/N?) and 0(1/Nf) respectively in this paper. 




(2.7) 



a = ^ + \r] 



(3 = l-7]-x 
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As a preliminary to introducing the formal equations which will be solved 
to give an expression for r\i we review briefly the leading order analysis. 
Whilst this was done initially in || that calculation was carried out in mo- 
mentum space. Although the same results will be obtained in coordinate 
space as in momentum space we will mostly work throughout this paper in 
coordinate space. The exponent r] is obtained by examining the skeleton 
Dyson equations with dressed propagators at the critical point g c , (2.3). As 
the fields obey asymptotic scaling there then one can replace the propaga- 
tors constituting the Dyson equations with (2.6) which will therefore result 
in equations involving a and (3 which can be solved. The Dyson equations 
which we consider here are illustrated in figs. 1 and 2 having been truncated 
at 0(1/N 2 ). (Ordinarily at 0(1/N 2 ) in this approach one has in addition 
several three loop graphs but these vanish in QED due to Furry's theorem.) 

The quantities and correspond to the inverse propagators and 
their asymptotic scaling forms are obtained from (2.6) by first transforming 
them to momentum space, using (2.7) and its derivatives, and then inverting 
them via G~ 1 G = 1, before transforming back to coordinate space, ffl, Qj. 



As in other critical point analyses of models involving a gauge field, [18], the 
asymptotic form of its inverse propagator is determined by restricting the 
inversion to be on the transverse subspace of momentum space, since only 
the transverse piece is physically relevant, [p|. Following this procedure, 



ij) 1 (x) ~ 

as x 2 — > where 

aa(a — fi) 



m((3) 
B(x 2 ) 2 »-P . 



r(a) 



TT 2fJ, (fi — a)a{a) 



, m(/3) 



2(2/x-/3) x fJi x u 
(20-2/x-l) x 2 

[4fc-/?) 2 -l]a(/3-/j) 
4vr 2 ^ - f3) 2 a(l3) 



(2.9) 
(2.10) 

(2.11) 



and (2.9) was first given in [0]. 

Thus with (2.6), (2.9) and (2.10) and, for the moment, retaining only the 
leading one loop graphs of figs. 1 and 2 the Dyson equations at criticality 
are equivalent to 



r(a - 1) + 



2(2/x-l)(/3-/i + l)z 
(2/i - 2/3 - 1) 



(2.12) 
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for the electron where we have factored off the common piece ^f. Also, 



_ ™>iP) 



(x 2 ) 2 ^~P 



2(2 fj, - (3) x^x u 

Vp> + (2/3 _ 2/U - 1) X 2 



4:zN f ( 2x^x u s 



( x 2)2«-l {Via* x 2 J 

(2-13) 

for the photon and we have set z = A 2 B. As the transverse part of the 
gauge field in momentum space is the only physically meaningful part of 
(2.13), [19], we project this out by first transforming to momentum space, 
multiply the resulting equation by the projection operator P fMU (k), before 
mapping back to coordinate space. Equivalently one can make the following 
replacement for the longitudinal components of (2.13) in x-space, valid for 
all exponents a, which corresponds to this operation, ie 

(x 2 ) a 2(a- l)^ 2 )"- 1 y ' 

Thus the relevant piece of (2.13) is 

(n-f3)m(f3) 4(q - l)zN f 
U " (2 M - 2/3 + 1) (2a -1) [ b) 

In writing down (2.15) we note that the powers of x 2 , which cancel in the 
leading order analysis of (2.12), are cancelled after projecting out the rele- 
vant piece of (2.13). 

We now have two equations, (2.12) and (2.15), involving two unknowns, 
z and rj, and so eliminating z between both the consistency equation for r/ 
emerges, ie 

(//-/?)(/? -/x + lM/?) , 2N f (a-l)r(a-l) 

[4(^-/3)2-1] + (2 M -l)(2a-l) [ - 0> 

which can be solved with a = [i + \r) and j3 = 1 to give 

( s (2 j »-l)(2-/z)r(2/,) 
m0l) = " 4T 2 ( P )T(p + in2-,) (2 ' 17) 

where r/ = J^i^i'Hi/^f an d hence from either (2.12) or (2.15), 

1 " 47rV(2/i-l)( / /-2) 1 8j 

which is required later. Although z\ = in three dimensions this does not 
imply that the formalism breaks down in this dimension. First, there are 



9 



0(l/N 2 ) corrections to z which ensure ^4^0 and B ^ and if one were to 
compute z\ in arbitrary covariant gauge then the factor becomes (2// — 3 + 6), 
[ |17f . Second, in the determination of the final results for exponents it turns 
out that the appearance of such (2/i — 3) factors always arise with a factor 
(2// — 3)~ x and therefore there is no difficulty in obtaining non-zero and 
non-singular results when the final d-dimensional results are restricted to 
d = 3. Indeed when (2.17) is evaluated in three dimensions it agrees with 
the wave function renormalization constant calculated in an explicit three 
dimensional large Nf renormalization of QED, f20(| , which is another strong 
check on our analysis in addition to the already stated consistency of the 
e-expansion of (2.17) with the leading order 0(g 3 ) terms of (2.14). Whilst 
(2.17) was obtained in a similar fashion in momentum space the coordinate 
space approach provides the starting point for computing the corrections to 
7/i. For instance, one can expand (2.16) to the next order in 1/Nt but this, 
of course, neglects the higher order graphs which have to be included. 



3 Corrections to consistency equations. 



In this section, we derive the formal corrections to the consistency equations 
to determine 7/2 by including the higher order corrections to the skeleton 
Dyson equations with dressed propagators. For QED there are only two 
such correction graphs which are given in figs. 1 and 2. First, we formally 
denote the values of the two loop integrals by £ for the electron self-energy 
and ILjj, = tj^H + x^x v E/x 2 for the photon self energy in coordinate space. 
We will discuss their explicit evaluation in subsequent sections but note that 
by value of the graph we mean the result one would obtain by computing the 
integrals with unit amplitudes and symmetry factors (such as minus signs 
for fermion loops etc) excluded. 

Thus, the graphs of figs. 1 and 2 are equivalent to 

= r{a - 1) + (x 2 ) x f{P)z + (x 2 ) 2 *z 2 £ (3.1) 



and 




m(J3) 
(x 2 ) 2 ^ 
z 2 N f 



2(2/j - (i) x^Xu 



(2/3 - 2/x - 1) 



4zN< 



f 



(x 2 ) 



2X nXp 



f x 2\4a+P-2ij,- 



(3.2) 



where f(/3) = 2(2//- l)(/?-/x + l)/(2/x-2/3- 1). Unlike at leading order we 
cannot cancel off the powers of x 2 since now x = 0(1/Nf) which will give 
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0(1/ N 2 ) contributions. Further, from the naive computation of £, II and 
H it turns out that they are in fact infinite due to divergences which arise 
from vertex subgraphs when x = 0. This situation, indeed, is completely 
analogous to that at 0(1/N 2 ) in other models, [|, ||. As a first step in 
treating these infinities we introduce a regulator, A, by shifting the exponent 
of the gauge field, (3, by (3 — > (3 — A, where A is an infinitesimal quantity 
playing much the same role as e = (4 — d) /2 does in dimensional regulariza- 
tion of four dimensional perturbative calculations. With the introduction of 
A we formally define the following quantities 



K P 

+ t! , n = 4- + n' 



(3.3) 



A A 

where the prime, denotes the completely finite parts of E, II and H with 
respect to A and both the residues, K, P and X, and finite pieces are purely 
functions of fi, a and (3. The poles with respect to A of (3.3) which will 
therefore appear in (3.1) and (3.2) are removed by the vertex counterterm 
which is available at each vertex in the one loop graphs of figs. 1 and 2. 
Thus denoting this counterterm by u, which can be expanded about u = 1 
in the large Nf expansion, as 



the regularized Dyson equations therefore become 


and 



r(a - 1) + zu 2 (x 2 )* +A f((3 - A) + z 2 (x 2 fx +2A ( ^ + £' 



(3.4) 



(3.5) 







m{[3 - A) 
4zu 2 N f 



2(2/1 -P + A) Xf^Xv 



{213-2/1-1- 2A) x 2 



2\2a-l 



2x ^x v 



z 2 N, 



f 



- P^ + X 



A K - lflu ' " x 2 



(3.6) 



With (3.4) and expanding each term of the electron equation, (3.5), to 
0{1/N 2 ) and the finite parts in A, then the divergent terms at 0(1/N 2 ) 
with respect to A are 

4u lZl (2n-l)(f3 - + 
(2(i -2(3-1) 



+ 



Kz 2 
A 



(3.7) 
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Setting (3.7) to zero gives a finite consistency equation for the electron so 
that the A — > limit can be achieved without difficulties. Of course this 
choice corresponds to a minimal scheme. If one were to absorb finite parts 
into ui then this would alter only the values of the amplitudes but not the 
exponents, [jl]]. The resulting finite equation, however, contains lnx 2 type 
terms which would otherwise spoil the analysis at criticality when x 2 — » 0, 
M. To avoid this difficulty, \i which has yet to be determined is defined in 
such a way that the lnx 2 terms are absent. So with 

the finite Dyson equation at 0(1/N?) in the critical region becomes 



r(a - 1) + 



f 

2z(2n - 1)(P - n + 1) 



(2/i - 2(5 - 1) 
Kz 2 

+ + (2 „- 2;3 _ 1)(/3 -, + 1) (3-9) 

Unlike in other 0(1/N 2 ) analyses, |2], |5| , where only the finite parts of 
the higher order corrections contributed to r]2, here, at least formally, the 
residue K appears in (3.9). This is due to the fact that in the non-regularized 
equation there is a non-zero function of (3 multiplying the term involving the 
counterterm, u. Thus when this function is expanded in powers of A, a finite 
term remains in the A — > limit, when the linear term of f(/3— A) multiplies 
the counterterm m, which involves the residue K. 

The treatment of the photon self energy corrections are somewhat similar 
to those of the electron, though as at leading order we consider only that 
part of (3.6) which corresponds to the transverse part in momentum space 
since this is the physically important piece. Thus, using (2.14) we consider 

2(fi - [3 + A)m(/3 - A) 8zu 2 {a - l)(x 2 )x+ A 







(2/j, - 2j3 + 1 + 2A)N f (2a - 1) 



2/ 2n,2 X +2A 



z (x 



n + 



(3.10) 



2(4a + P - 2/i - 2 - A) _ 
where we note that the same counterterm as (3.5), of course, arises. Again 
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isolating the divergent terms with respect to A at 0(1/ N?) they are ab- 



sorbed by setting 

(2a - l)zt 

U\ = — : 

16(a - 1) 



2(4a + (3-2^-2) 



(3.11) 
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which does not appear to be equivalent to that obtained from (3.7). For 
the moment, we note that when the explicit values for P, X and K are 
determined in a later section, we will observe that (3.7) and (3.11) are in 
agreement. Further, In a; 2 terms are removed by again defining xi appropri- 
ately. In this case, we set 



Xi 



(2a - l)zi 
8(a - 1) 



P + 



X 



2(4a + p-2p-2) 



(3.12) 



which will, of course, also agree with (3.8) giving us at least one check on the 
explicit evaluation of E and 11^,. Consequently, the finite Dyson equation 
for the photon at criticality is 







2(ji - P)m(P) 8z(a - 1) 



— z 



(2fi -2P + l)N f 

n' + 



(2a 



+ 



-1) 
X 



2(2a-l) 2(2a-l) 2 . 



(3.13) 



where, like (3.9), the residue of the two loop correction, X, will also give a 
contribution. It arises in a different way to the appearance of K in (3.9), 
through the expansion in powers of A of the coefficient which multiplies S 
after restriction to the transverse piece. 

Whilst we have formally derived the finite corrections to (2.12) and (2.15) 
we require the explicit values of E, II and S, and hence xii which is needed 
to evaluate the functions of P in (3.9) and (3.13) to 0(1/Nj) in order to 
obtain the formal consistency equation satisfied by 772- 



4 Computational tools for computing £ and li^. 

Before detailing the explicit calculation of the two loop corrections we will 
now review and develop the necessary techniques which will be required. 
First, we recall that in solving the bosonic and supersymmetric 0(N) a 
models at 0(1/ N 2 ) in Jl], 0, |5|, M, extensive use was made of the technique 
known as uniqueness which was first introduced in |2lJ] and subsequently 
used and extended in various forms in [22-24]. It is applicable to models 
which involve a 3-vertex, where the exponents of the propagators forming 
the vertex are initially arbitrary. In endeavouring to compute such a 3-vertex 
in coordinate space it turns out that the calculation cannot be completed 
in closed form unless the sum of the exponents are restricted to be the 
dimension of spacetime, 2fi, for a purely bosonic vertex with no derivative 
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couplings, pi] ]. When this uniqueness criterion is satisfied then the integral 
can be completed to yield the product of various propagators multiplied 
by a factor dependent on the exponents of the initial vertex. We have 
summarized this rule in fig. 3 where the Greek letter beside each line denotes 
the exponent of that propagator and we have written the product of the 
three resulting propagators graphically as a triangle. Also, v[a\, a.2, 03) = 
7r M n?=i a ( a i)- Further, if a + j3 + 7 = 2/i, + n for any positive integer n 
then this vertex is also integrable, yielding a product of triangle graphs. We 
have outlined this rule in detail since in the a models its 3- vertex exponents 
indeed sum to the uniqueness value, [Q, ||], which meant that the technique 
could be applied there. More recently, a similar rule was developed for the 
vertex of the Gross Neveu model, where the basic uniqueness condition due 
to the presence of fermions becomes 2\i + 1, |lfi]. 

In the model we are concerned with here the relevant 3-vertex contains 
a gauge field interacting with two fermionic fields and therefore we need 
to develop the analogous uniqueness integration value for this vertex. It is 
illustrated graphically in fig. 4 in coordinate space, where the indices \i and 
v refer to those which appear in (2.6) and the numerator of the integral is 
{i ~ ~ ?0- We use the convention that a fermion propagating from 

x to y has a factor — f). Repeating the analysis described to derive the 
result of fig. 3 yields a more involved expression. After introducing Feynman 
parameters for each propagator and completing two of the four integrations 
one obtains a sum of various integrals each involving one hypergeometric 
function. Essentially, the uniqueness condition emerges by choosing the 
individual arguments of this function in such a way that it is equivalent to 
a simple algebraic function, after which the integral can be computed and 
the result of fig. 3 obtained for the purely bosonic vertex. In the case of 
fig. 4, it turns out that the minimal uniqueness condition for this vertex is 
a i + Oi2 + «3 = + 2, for the Landau gauge, which results in a sum of 
triangle graphs analogous to the right side of fig. 3 but with bosonic and 
fermionic propagators. We do not give the full graphical expression here 
for the following simple reason. If one now examines the electron photon 
vertex, then the value of the vertex, when the leading order exponents are 
substituted, is 2a + f3 = 2(j, + 1, which is one step from the uniqueness 
value unfortunately. This means that unlike the models of [||, ||], g], we do 
not have a direct integration rule with which to perform our calculations. 

To circumvent this difficulty a different approach will be required which 
we now outline. The ideas we use, however, derive from similar difficulties in 
other models. For instance, in the Gross Neveu model, Q , it was noted that 
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one can always rewrite a graph involving fermions in terms of its constituent 
purely bosonic graphs, by explicitly computing the trace over 7-matrices. A 
consequence of this is that some of the exponents of the resulting graphs 
are reduced by an integer so that the vertices become equal to the (bosonic) 
uniqueness value and therefore integrable from fig. 3. Whilst it was possible 
to introduce methods in the Gross Neveu model, where taking the trace 
is not necessary, it is an important observation for the present case. The 
application of this to QED will, in principle, be similar though to derive the 
constituent bosonic graphs involves taking traces over a larger number of 
7-matrices. We will mention further simplifying aspects later but note that 
in our manipulations we will make use only of the algebra (A.l) which is, of 
course, valid in arbitrary dimensions. 

Secondly, one technique which is used frequently in computing massless 
Feynman integrals is that of integration by parts, [0]. Indeed a recursion 
relation was developed for a two loop integral in || which was primarily 
required for changing the value of the vertices by ±1, to obtain two loop 
graphs which had either integrable triangles or vertices. Equally we can 
apply this technique to the gauge vertex of fig. 4 where the gauge field 
propagator is 

1 

and its numerator was noted earlier. In the corresponding integral we choose 
to rewrite part of the second term of (4.1) as 

z fi 19 1 

so that when one integrates by parts, factors from differentiating the denom- 
inators will combine with the numerator to reduce the number of 7-matrices 
present by two in some terms. Thus carrying this out and rearranging, one 
obtains the result illustrated graphically in fig. 5 for arbitrary values of the 
exponents. Therefore, one can rewrite the 3-vertex involving a photon in 
Landau gauge in terms of three other graphs which have a simpler structure. 
For instance, the first graph is proportional to what one would obtain if the 
second term of (4.1) was ignored. In the remaining two 3- vertices the long 
dashed vertical line, with an index at its external end, joining the point of 
integration, z, is to be understood as the propagator z fl /(z 2 ) ai , and we also 
note they involve only one 7-matrix. The usefulness of this rule will become 
apparent later but we remark that using it in a two loop graph will, for 



2e*i 



(2u-2qi - 1) 



(4.1) 
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instance, alter the structure of the propagator joining to the external ends 
of the original 3- vertex, in the case of the second and third terms of the right 
side of fig. 5. Therefore, two further rules are required, which are illustrated 
in figs. 6 and 7, and each is established by an integration by parts similar 
to (4.2). We note that the numerator of the graph of the left side of fig. 6 
is (J) — — whilst that of fig. 7 is {jj — 9^)7^. 

5 Computation of II^j,. 

In this section, we outline the calculation of Ii lJil y(x) and concentrate on 
the main features. As discussed previously, since the method of uniqueness 
cannot be used directly, we have to develop an algorithm which will give rise 
to integrals to which one can apply this method. 

First, the detailed graphical expression for Ii llll {x) is given in fig. 8, 
where the Greek indices at the vertices (internal and external) correspond 
to the Lorentz index of the corresponding uncontracted 7-matrix. Therefore, 
this two loop integral involves a trace over eight 7-matrices. As a first step 
in reducing this number, we use the integration by parts rules given in the 
previous section. Applying fig. 5 first yields three graphs, one of which 
is the first graph of fig. 9, with the remainder treated using fig. 6 at the 
opposite vertex. After a suitable rearrangement they correspond to the final 
three graphs of fig. 9 where one has, of course, to include the regularization 
A. We have given only the graphical form of these graphs since the simple 
factors associated with each can be readily deduced from figs. 5 and 6. The 
final graph is multiplied by the factor, 

\x{$f£f) = 4[2xV - x 2 rf v \ (5.1) 

For simplicity, we denote the contribution from the first graph of fig. 9, with 
correct factors by / and that from the others by //, treating each separately. 

To illustrate the type of calculations required once the integral has 
been broken up into its constituents, we concentrate for the moment on 
the final graph of fig. 9, ie (a, a, a, a, (3 — 1 — A), where we have defined 
(ai, «2 ; «3 ; «4> «s) to be the general bosonic two loop integral of fig. 10. To 
compute (a, a, a, a, (3 — 1 — A), we use the methods of subtractions devel- 
oped in Q and note that in the absence of the regulator, A, the graph is 
integrable using the uniqueness rule of fig. 3, since 2a + ((3 — 1) = 2/x, 
as a = fi and (3 = 1, at this approximation in large Nf. Carrying out the 
naive integration at A = 0, will yield a result proportional to T(0), which 
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is infinite. However, when the theory is regularized by A, one loses unique- 
ness at each vertex and cannot naively apply fig. 3. Instead, recalling that 
only the pole and finite pieces of IT^,, with respect to A, are required to 
deduce 772, one need only consider the combination {A — B) + B, where A = 
(a, a, a, a, (3 — 1 — A), and B is such that in the presence of A it is calculable 
but has the same divergence structure as A. Thus, (A — B) is finite and can 
be deduced for A = 0, |^. A suitable choice for B is 

(a,a,a,0,/3 - 1 - A) + (a, a, 0, a, f3 - 1 - A) (5.2) 

which corresponds to a sequence of integrable chains of propagators for A 
7^ 0. The combination (^4 — B) is obtained by completing one integration 
with A = first, before shifting each exponent of the propagators of the 
resulting chain by a temporary regulator, 5, so that poles in 1/5 can be seen 
to cancel before setting 5 to zero, J3]. Thus, one obtains 

a 1 m 2A 2 (a)a(/3- 1) 
{a, a, a, a, (3 — 1 — A) = rV/AA 



1 ((3-1) 

+ TT, 7T + 



1 + A (B(f3) - B(a — 1) 
(5.3) 



(0 -1) 2(a-l)( M -/3) 

where B(x) = ipdi — x) + vp(x), and ^(x) is the logarithmic derivative of the 
T-function, without specifying the leading order values of a and (3 as fi and 
1, respectively. This last remark is actually crucial, since naively setting a 
= fi and (3 = 1 in (5.3) one immediately discovers that the formal expression 
appears to be infinite which is due to the fact that the value a = fi is the 
anti-uniqueness value of a bosonic field, f2|. In reality this problem is not 
significant since there are cancelling infinities which arise from other graphs 
comprising 11^, so that the sum of all the constituent pieces is indeed finite 
upon setting a = fi and (3 = 1 in the final expression. Indeed this cancellation 
will provide us with a stringent check on our analysis and was also a feature 
of the much simpler Gross Neveu model, j|, if one rewrites the 0(1/N 2 ) 
corrections there in terms of its bosonic components. 

For the remainder of the calculation one essentially breaks up the origi- 
nal integral into similar bosonic two loop graphs and computes them using 
this method of subtractions though some are not as straightforward as the 
one outlined above. For instance, the first two graphs of fig. 11 can be 
treated directly, except that in the second there is only one subtraction to 
consider, since the divergence arises in the left vertex subgraph. Whilst 
the vertices of the third graph of fig. 11 are unique, one cannot apply the 
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subtraction procedure directly, since the infinity structure of the naive sub- 
traction (a, 0, a — 1, a, j3 — A) does not match that of the original. Instead 
an alternative technique is required which involves rewriting the integral in 
terms of other graphs, which can be computed by subtractions. In [Q], recur- 
sion relations for graphs of the type of fig. 10 were developed by integration 
by parts which have the effect that the exponents of the lines comprising the 
new integrals are adjusted by ±1 relative to the original. Whilst we were 
unable to make use of that rule here it is clear that a similar rule will be 
needed. 

Such an alternative was given in [24, pH], being derived in the latter by 



considering the uniqueness rule for a bosonic 3-vertex when its value is 1\x 
+ 1 and is illustrated in fig. 12. Applying it to the upper internal vertex of 



the graph of fig. 10, for example, one obtains the recursion relation, [25], 



/ \ »50-«5-l) , -t -| , 1V 

a 2 , a 3 , a 4 , a 5 ) = — - — — ,(ai - 1, a 2 - 1, a 3 , a 4 , a 5 + 1) 

(«i - l)(a 2 - 1)^ 

(ai+ag-M-l) , -. . 

+ ; 7T (ai - 1, 02,03,04, «5/ 

(ai-1) 

(ai + «2 — /i — 1) , 1 , 

+ -, rr (ai,a 2 - l,a3,a 4 ,a 5 ) (5.4) 

[p.2 ~ 1) 

Further, one can derive more relations by either applying fig. 12 to other 
vertices or by first making use of the possible general transformations given 
in the table of [Q], then applying the rule before undoing the initial trans- 
formation by applying its inverse. In appendix B, we have listed the other 
recursion relations we required for this calculation, but note that we believe 
this list is not exhaustive. Returning to the third graph of fig. 11 and ap- 
plying (B.10), results in the following three integrals, (a, a, a, a, (3 — 1 — A), 
(a, a — 1, a, a — 1, (3 — A) and (a — 1, a, a — 1, a, (3 — A), where the final two 
are equivalent under a change of integration variable and each is calculable 
by the subtraction procedure. 

In the first part of appendix C, we have listed the basic library of two loop 
building block integrals we required to compute 11^^, each expanded to its 
finite part with respect to A. However, we have not given the values of those 
integrals which can be derived directly by subtractions, only those which 
used the recursion relations given in appendix B since they are more tedious 
to deduce. Where possible, we have calculated several of these integrals in 
independent ways, using different recursion relations and so are confident 
that the expressions are in fact correct. 
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One basic integral, (a, a, a, a, /3 — 2 — A), deserves special attention, 
since it does not have any uniqueness values at A = and also cannot be 
related to any other basic integral for reasons which will become apparent. 
We compute it by first applying the transformation <— , in the notation of 
[0, and then use the recursion relation (B.12), which results in (a,/i — a + 
1, fi — a, a, fj, — 1 — A) and {a — 1, fj, — a + 1, — a, a, fj, — 1 — A). Applying <— 
to the former yields (a, a, a — 1, a, j3 — 1 — A) of which only the residue at 
the pole with respect to A is relevant since the associated factor contains A. 
One then applies (B.8) appropriately to the other graph before undoing the 
initial <— transformation. One of the resulting two graphs can be treated as 
in an earlier discussion whilst the other is (a — 1, a — 1, a — 1, a — 1, j3 — A). 
We were unable to evaluate this graph directly for arbitrary a and /3, though 
it is in fact finite, both with respect to A and setting a = fi and = 1, 
and had therefore to leave it unevaluated in our analysis until we set a = 
\i and (3 = 1 after all contributions to had been summed. After this 
substitution we can relate it to a known integral ChT(l, 1), in the notation 
of ||, by applying the transformation — >. Thus, 

(/i-l, M -l,M-l,M-M) = "(^ Z 3) ChT ^ V ( 5 - 5 ) 

and ChT(l, 1) = 3vr 2 ^a(2/i - 2)T(^ - l)[ip'(n - 1) - ^'(1)], @. Adding all 
the pieces, we have 

2^(^-l)V(a)a(/3-l) 
(a,a,a,a,/3-2-A) = (a _ 1)2(j g _ 1)rM 

2Q3 + 1 -//)(/? -/i)(2a- 3) 
(a - l) 3 

x (a — 1, a — 1, a — 1, a — 1, /3 — A) (5.6) 



+ 



Several other basic building block integrals also involved (a — l,a — l,a — 
1, a — 1,0 — A), which have therefore been left unevaluated in the list in 
appendix C. 

We conclude the section by returning to the graphs of fig. 9 and briefly 
mention the key points required to compute the first three aside from those 
already discussed. For the remaining graphs of II, one occurs in the evalu- 
ation of /, whilst the other can be treated by a subtraction, even though it 
involves one fermion propagator. (The rules for integrating chains of prop- 
agators involving fermions were given in 0, for example.) 
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Finally, to treat / we have used the result 

- 27V7 ff - 2(/i-l) 7 Y/ (5.7) 

which is valid for arbitrary dimensions and has been derived solely from 
(A.l). It turns out that the three terms of (5.7) with three 7-matrices give 
graphs, which although involve fermions, can be computed by the subtrac- 
tion procedure without having to explicitly project out the Tf" v and x^x v 
components. For instance, the graph involving only the three 7-matrices of 
the last term of (5.7) is 

8A 2 (a - l)a(f3)Q^ 



A(a-l)y(/i-l) 



1 + A B(J3) - B(a - 1) - 



1 



(5.8) 



(a-1) 

where Q^uix) = f]^ — 2x^x u /x 2 . For the graph corresponding to the first 
term of (5.7), we had to treat its transverse and longitudinal components 
separately but in doing so each piece involved at most four 7-matrices, and 
therefore we employed (A. 5). 

Thus having given an extensive discussion of how to compute Ii^ v we 
note that adding all contributions, i" is, at a = n and (3 = 1, 



32tt 2 ^(// 
Ar2( M )(2 M -3) 




2 36(/i) 
//(2/i-3) ~ 2(7* - 1) 



where <->(//) = - 1) - ^(1). Whilst we find 



2A 



(2/x - 3) 



(2/x-3)r2( M )A 
for II. Thus summing (5.8) and (5.9) we obtain 

16vr 2 ^ ( 2 Xfl x v 



(5.9) 



(5.10) 



n 



flU 



+ 



Ar 2 ( /U )(2 / u-3) 
Af3(/x-l)6(/x) 



(2/x-l)(2- / u) 



+ 



2(2(1 



As mentioned earlier we have used the stringent check that the result for 
must not be singular at a = (i and (3 = 1 when adding all the contributions 
from all the pieces. Also the tensor structure is the same as at one loop and 
if one were to transform (5.11) in the context of (3.6) to momentum space 
the transverse projector, P^ u (k), will emerge so that the result is gauge 
invariant. This provides another useful check on the calculation since the 
individual constituent graphs do not each have the necessary Q^ v structure. 
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6 Computation of E. 



In the previous section, we outlined in detail the methods required to de- 
termine H^ u . The algorithm discussed there can also be used to calculate 
E though there are several differences which we mention here. First, in E 
of fig. 1 we have to compute a graph involving two gauge fields. Again the 
first step is to break the graph up by using the integration by parts rule on 
one of the gauge fields which will result in three graphs, denoted by /, II 
and III, respectively corresponding to the three terms of fig. 5. To proceed 
further one applies either the rule of fig. 5 again or that of fig. 7. 

It turns out that the resulting nine graphs are much more straightforward 
to analyse in terms of the basic bosonic two loop building block integrals than 
for Ii^ u and in particular there is no graph analogous to (a, a, a, a, f3— 2— A). 
However, as the graph E will be proportional to $ overall, where we use the 
convention that the left external vertex of fig. 1 is the origin, the result 
of any manipulations will yield graphs with at most three 7-matrices using 
the results (A.2)-(A.4). To proceed further, one can either, if possible, use 
fermion subtractions similar to the bosonic subtractions of section 5, which 
were introduced in Q, or, if, for example, the naive subtraction is insufficient 
to match the divergence structure of the particular graph then it is best 
to take the explicit trace by first multiplying the graph by ^. This will 
yield a sum of purely bosonic graphs which can be treated individually. We 
have given the library of essential integrals in the latter half of appendix C 
where we have again listed only those expressions which were derived from 
recursion relations. The two loop graphs which are directly obtained by the 
same type of subtractions discussed in the previous section are not given 
since they are easy to deduce. 

As an example we discuss the computation of the graph obtained after 
applying the integration by parts rule of fig. 5 to each of the internal vertices, 
which is equivalent to that of fig. 1, but where the second term of (4.1) is 
absent. For concreteness, the numerator of this integral is 

l v {-f)l°{i-f)lAt-t)l« (6.1) 

where y is the location of the top vertex of integration and z the lower which 
involves seven 7-matrices. However, this number can be reduced by using 
(A. 4) several times and gives after a suitable rearrangement 

4(v-l)M(x-y) 2 -(y-z) 2 -(x-z) 2 } + (/i - 4)M - j)U - (6.2) 
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where we have made a change of variables to obtain the factor 2 in the 
first term. Whilst the final term involves three fermion propagators one 
observes that its form is equivalent to the fermion self energy of the Gross 
Neveu model, which was calculated in Q. If we introduce the notation 
(di, «2> "3; «4) "5) to mean that the propagator with exponent o>\ corre- 
sponds to a fermion then (a, (3 — A, a, (3 — A, a — 1) can be calculated directly 
by using the subtraction (a, (3 — A, a, 0, a — 1) and the result, to the finite 
part in A, is 

n 2 »a 2 (a - l)a{(3) 



(6.3) 



AT(n)(n -a)(a- 1) 

Again, one observes that this expression is infinite at a = n and (3 = 1 but 
in the overall sum for S, such infinities will in fact cancel. For the remaining 
graphs, (a, /3 — A, a — 1, f3 — A, a) and (a, (3 — 1 — A, a, (3 — A, a), one has 
to take a trace for each to obtain 



27r 2 ^a(a - l)o(a)a(/9) 



AT(/z) 

for the former and 

7r 2/i a(a)o(a — l)a(/3) 



l + A(B( / 9)-B(a-l)-i+ (/3 " 



^ 4/x(a — 1) 



(6.4) 



+ A 



2Ar(/i) 

(0 



6 



(2 M + i)(/?-i; 

/i(a — 1) 



,M W )-B(a.-1)] + ^^ 

(a — 1)^ 



//(a — 1) 

(/3-l)(yu-l) (/?-l)(^ 2 -3/i + l) 



4 
/i 



/i(a — 1) 



^ 2 (a-l) 



(6.5) 



for the latter. The computation of the remaining graphs is also similar. 

Finally, we close this section by giving the results of the three component 
graphs which comprise S, as an aid to the interested reader, having checked 
that each is finite after the substitution of a = /i and (3 = 1. Thus 
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16tt 2 ^(2^ - l)(/i - 1) 

A(2/x - 3)Vr 2 ( M ) 

87^(2/1 - 1) 
"A(2/,-3) 2 r 2 (^) 
4tt 2 ^(2^ 
(2/,-3)Vr 2 (/,) 



r V 2/Lt (2/x-3) 

4(M-5 



/i-2 + A - 



(2/x - 3) 



(6.6) 
(6.7) 
(6.8) 
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We note that one of the constituent graphs of T,jj involves the relatively 
large expression (C.17) which is deduced from a recursion relation containing 
(C.15) and (C.16). Hence, 



+ A (2/i - 5)// + 



2(2/x-l)( /J -2) 2 

(6.9) 



/jr 2 (/j)(2/j-3) 2 A 

4(/x-l) 2 ( M -2) 8(2 / u-l)(/,-2) 2 



// (2/x - 3) 



We conclude our discussion of the two loop calculations by remarking 
that in computing both and £ we made extensive use of the relation 2a 
+ /3 = 2fj, + 1 between the exponents a and /?, which is valid at the order 
we are calculating, to simplify substantial amounts of tedious algebra. 



7 Derivation of 772. 

Having discussed the derivation of the two loop corrections we now return to 
the formalism developed in section 3 and derive the consistency equation for 
7/2. From (5.11) and (6.9), we first of all check that the values one obtains 
for xi m both renormalizations agree. Thus with 

87^(2/, -l) 2 (^-2) 2 167^(2/1 -l)(2-/z) 

(2/, - 3)2r( M )rv + 1) ' (2/x - 3)r( M )i> + 1) 1 ' J 

then 

Xi = - m (7-2) 



from either (3.8) or (3.12), and, as was noted in |26J, this corresponds to the 
QED Ward identity. Likewise, the expressions for the vertex renormalization 
constant, u, both agree. With this value for xi, then (3 = 1 + 0(1/N^), so 
that we can now write down the formal consistency equation for 772, which 
is given by eliminating z<i from (3.9) and (3.13), as 

2r?2 1 1 

r,\ /x (2/i-l)( M -l) 

(2/x-3) 2 r(/i)r(/x + l) / K 



4tt^(2/x- l) 2 (/j-2) 2 V (2/x-3)(/i-2) 
(2/x - 3)iW( M + 1) / 5- X \ 

167t^(m-1)(/x-2) v 2(2/x- 1) 2(2/x- 1)V 1 ' ' 
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The explicit expressions for II' and £' can be read off from (5.11) and (6.9), 
respectively, and with (7.1) some straightforward algebra therefore leads to 



2 [ 3//(/i - l)9(/x) , 3 , 3 



(2^-l)( M -2) 2/x (m-1) 
1 28 35 



(7.4) 



3(/x - 2) 2 9(2^ - 1) 18(/i - 2) 



which is the main result of this paper, and is an arbitrary dimension ex- 
pression for the 0(l/Nf) part of the elect ron anomalous dimension in the 
Landau gauge and as was noted earlier encodes information on the gauge 
independent part of the 4-dimensional renormalization group function, "f(g). 

Aside from the internal consistency checks on the explicit derivation of 
the values of the two loop integrals already discussed, we have checked that 
the final result for 772 agrees with the three loop MS anomalous dimension, 
j(g c ), expanded near four dimensions. It is worth noting that this three loop 
result is in fact another very stringent check on our result since we have only 
evaluated two 2-loop graphs. Moreover, we can now derive the coefficients 
of the higher order terms of 7(5) at 0(1/Nj), which have not been given 
previously. For example, from (2.5), with (2.3) and (7.4), we find 



where the unknown constants c\ and d\ can only be deduced from 77 at 
0(1/Nj) and 0(1/Njr) respectively. The next to leading order coefficients 
with respect to Nf, at each subsequent perturbative order, are deduced by 
first extracting all the 0(1 /Nf) corrections to g c , which are contained within 
the large Nf (3- function of ||. Comparing the numerical structure of (7.5) 
with the four loop /3-function of (2.2), we note that they are similar in that 
the transcendental number, C(3), appears at fourth order but not at leading 
order in Nf. 

As (7.4) is valid in arbitrary dimensions, we can evaluate it in three 
dimensions and find 



7(5) 



+ 




+ c 1 N f + d 1 g 4 + 0(g 5 ) (7.5) 



rj = - 



3ir 2 N f 



8 



+ 



16(32 - 3vr 2 ) 
9vr 4 A^ 




(7.6) 
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where the 0(1/Nf) term agrees with the calculation of p0| ]. From (7.6), we 
can obtain estimates for the anomalous dimension for various values of Nf. 
For instance, at Nf = 3, r/ = — 0.09 and at Nf = 4, rj = — 0.06. 



8 Discussion. 

We conclude with various observations. First, we have given an analytic 
expression for the electron anomalous dimension at 0(1/N^) by solving the 
appropriate Dyson equations at criticality within the large Nf expansion, 
algebraically. Indeed, the information contained in the d-dimensional ex- 
pression, (7.4), relates to the perturbative coefficients of the gauge inde- 
pendent part of the electron anomalous dimension. For completeness, we 
note that this model is now solved at leading order in large Nf, since the 
critical exponents, or equivalently, the appropriate renormalization group 
functions, relating to the /^-function and the electron mass anomalous di- 
mension, 7m(<7), have been given elsewhere, ||, and in the notation of this 
paper, they are 

7 - fc) = - (p - Vf <fU) 

Nf 

where 2A = —f3'(g c ). So, for example, in three dimensions there is no 
0{\/Nf) correction to the /3-function which is gauge independent and this 
is consistent with the unltraviolet superrenormalizability of that theory. In- 
deed the four dimensional perturbative information, which is encoded within 
(7.4), (8.1) and (8.2), will provide useful checks for future explicit pertur- 
bative calculations. Second, with the techniques given in this paper, it 
ought now to be possible to go beyond the 0{1/Nf) expressions in (8.1) and 
(8.2), by extending our critical Dyson equation approach. For instance, the 
0(1/N 2 ) corrections to the /3-function of the bosonic and supersymmetric 
O(N) a models are known, [||, ||]. Indeed, the method of |l], || is such that 
one has always to first compute the exponent 77, which we have done here, 
before attempting to determine A at the same order in large N since, for ex- 
ample, one needs to know Z2 which is determined from (3.9) and (3.13) once 
7)2 is available. Finally, we remark that the methods which have been devel- 
oped here for computing anomalous dimensions in an abelian gauge theory, 



25 



will be applicable in solving similar models, such as the bosonic CP(N) a 
model, beyond the 0(1/N) exponents which are presently known, p7| . 

Acknowledgements. The author thanks Dr D.J. Broadhurst for providing 
a copy of reference |l5| and for discussions on the gauge independence of 
7(5) beyond one loop in QED, as well as Dr O.V. Tarasov for an indirect 
communication giving the three loop structure of 7(5). Drs D.R.T. Jones 
and H. Osborn are also thanked for brief conversations. 



26 



A Conventions. 



In this appendix, we briefly list our conventions and several results involving 
7-matrices which we required. Although we work in arbitrary dimensions 
we take the trace convention to be trl = 4, and to manipulate 7-matrices 
we made use only of 

{ 7 ^} = 2rT (A.l) 
working in Euclidean space throughout. From (A.l) we derive, 

7 CT 7 M 7a = -2(/*-l)7" (A.2) 
f7"/7, = V" + 2(/i- 2) 7 'V (A.3) 
y^y/y^ = - 27^7 CT 7^ 

- 27" W - 2(/i - lJWV 7 (A.4) 

Although we were able to manipulate the two loop integrals to reduce the 
number of 7-matrices involved, we needed 

tr(yyW) = 4[r^V P - rfrf" + rf p r} V(J \ (A.5) 
to complete the calculations of ELj, and E. 



B Summary of recursion relations. 

To compute various basic bosonic two loop building block integrals making 
up £ and 11^, we required various recursion relations, which we list in this 
appendix. Whilst this list is perhaps not complete, given the large number 
of tranformations one can make on the basic two loop integral of fig. 10, they 
were all that we required for our purposes. To compactify our expressions 
a little, we define a new five argument quantity similar to (01,02,03,04,05), 
but involving square brackets. In the following it will denote the two loop 
graph obtained from the general graph of fig. 10, but with the respective 
exponents adjusted by the arguments of the expression involving square 
brackets. For instance, [0,0,0,0,1] means {aii, 012,(23, ($4,015+1) etc. Also, 
when a square bracket appears with a subscript ± it will correspond to 
multiplying the overall expression by a factor (x 2 )^ , so that each term has 
the correct dimensions. We also use the notation of |2j and define 

0\ + 2 + 05 = Si , 03 + «4 + 05 = S2 
0\ + O4 + 05 = t\ , 02 + O3 + 05 = t2 

o\ + 02 + 03 + a 4 + a 5 = d (B.l) 
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From the rule of fig. 12, we have derived the following recursion relations, 
which we believe have not been given before, and emphasise that the rela- 
tions are completely general and applicable to any graph of the form of fig. 
10, and not solely for the specific problem dealt with in this paper. So, we 

have, 

(3/i-d)(d-2/i-l). 
(0:1,02,0:3, 04,0:5) = — - - - [0,-1,-1,0, 1J + 

+ ^±^^[0,-1,0,0,1] 

+ ^±^^ [0,0,-1,0,1] (B.2) 
(03 - 1) 

~ ^^11.-1.1.0,-1] 

+ (a2 ^ 5 _,_ 1) [0l _ MM1 

+ ^±^^ [0,0,1,0,-1] (B.3) 

(«5 - 1) 

- ^(M-«5-l) r -1 j _i j0 ,0,1] + 



(01 - l)(o 2 - 1) 



+ (ni ^-T' ) [0 --'- -°- 01+ 
+ (Q1 + g2 p-l) 

(01 - 1) 

5 (/X-0 2 -Q 5 -l) _ 

(2^- Sl -l)( Sl -/,) [U ' U ' U ' MJ 

02(^-02-05-1) 

+ (2 M - Sl -l)( Sl - /i ) LU ' 1 ' U ' U ' UJ - 

+ 7^ ^7 ^[-1,1, -1,0,1] (B.5) 

(2^-si - l)(si -n) 

Os(2m-*i)(*i-M-1) .[_i, ,1,0,0] 



(2/u - s 2 - 1)(«2 - At)(oi - 1) 



05(01- 03 -l)(M-«5-l) f _ , 

+ o 3 (oi-o 3 -l) [0 1 0] _ (R6) 
(s 2 - //)(2/x - s 2 - 1) 
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«3(*2-M)(M-«3-l) [0)1)1)0j -1]. 



(d-2jLt)(3/i-d-l)(a 5 -l)' 



, - a 2 - Q=3 - IKM ~ «4 ~ 1) r i n i _ii 

+ (d-2/z)(3^-d-l)(a 5 -l) L J ~ 

[0,1,1,0,-1]- 



(2/x-ti)(t 1 - M- 1) 



a 3 ( M -a 2 -a3-l) 
+ (2 M -t 1 )(t a - M -l) L ' ' ' ' J 

(2/i - - /i - 1) 
(2 M -^)(2 M -^) [0)0)0)05 _ 1]+ 



(2 M -t 2 )(t 2 - M -l) 



+ 0g aa )( d + aB - 3M -l) [0 Q 1>0 0] 
(2/i - i 2 )(*2 — M — 1) 
(2^-si)(d + a 5 -3AX-l) rn , , . 

+ ( 2/1 -t 2 )(t 2 -^-l) t '" 1 ' ' ' ^ ( ) 

(2 M -ti-l)(ti-^) 
(« 2 -l)(a 3 -l) t '- 1 '- 1 ' ' 1 ^ 

+ (^^pl) [0 ,-i,o,-i,i ]+ 

+ (a2 +;;; 1 / ;" 1) [-i,o,-i,o,i] + ( b.i 0) 

a8(t a -M)(2/i-l-<S,) r. 1>0>1)0j0] 



(ai - l)(*i -/u- l)(2//-t 1 ) 1 



(qi - a 3 - l)(3/i - ri)(rf - 2/z - 1) 
(ai - l)(ti - n - l)(2/i - ii) 

- J^^T^ A-A (B.11, 

Further, we record the recursion relation derived using integration by parts 
which we also required in evaluating H^ u , which is, g, |22| , 

= ^([0,1, 0,0,0]- -[-1,1, 0,0,0]) 

+ asQO, 0,1,0,0]_ -[0,0, 1,-1,0]) (B.12) 
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C Basic two loop integrals. 



In this appendix, we give a list of the basic bosonic two loop integrals re- 
quired to compute E, II and S, expanded to the finite part with respect to 
A. As discussed in sect. 5 these fall into classes ie those which are computed 
directly using the method of subtractions and those which are not but which 
are determined by using recursion relations. As the former set are easy to 
establish, we list only those of the second class which we require. First, we 
consider the basic graphs for II and S, using the notation of fig. 10. We note 
that the integral (a, a, a, a, [3 — 2 — A) has already been discussed earlier. 



(a, a — 1, a — 1, a, f3 — A) 



2tt 2 ^(/h - l)o(a)a(a - l)a((3) 

rV + i) 



(C.l) 



(a 



1, a, a, a, (3 — A) 



2^a{a)a{a - l)a((3) 



fl(a-l) 



^ + B((3) 



((3 - 1) 2/i(a - 1) 



(C.2) 



(a — 1, a — 1, a, a, (3 — A) = 



27T 2 ^a(a - l)a(a)a(/?) 



1 2 

- Ba-1 - + T 

a — 1 M — 1 



(C.3) 



(a, a — 1, a — 1, a, (3 — 1 — A) 



7r 2 ^a(/3 - 1) 
a 2 (/i-a)r(//)|A 



+ B(j3)-B(a-1) 



+ 



+ 



(/? + l-/i)(2a-3) 
(a - l) 2 



x (a - l,a - l,a - l,a - (C.4) 



(a — 1, a, a — 1, a, (3 — 1 — A) 



+ 



tt 2 ^- l)a 2 (a)a(/?) 
(/x-l-/9)(2a-3) 



(a - l) 2 

x (a — 1, a — 1, a — 1, a — 1, (3) (C.5) 
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(a, a, a — 2, a, (3 — A) 



Tr 2 »a 2 (a)a(P - 1) 



+ 
+ 



rv + i) 

H 2 + n - 1 3 (g - l) 2 
M/i-l) ~ 2 ~ (//-/?) 
a a(a- 2)(/3 — 1) 



I + 5(^)-S(a-l) 



2(a- l)(n-0) 



(C.6) 



(a, a — 1, a — 2, a, (3 — A) = 



tt 2 ^(2 - n)a(a)a(a - l)a(/3) 



r(/* + i) 

H 2 + fi-l fi-1 



- B(a-1) + 

09 -1) 
(/z-2)(a-l). 

Next, the basic integrals for £ are, 

ir 2 v a (a - \)a(a)a((3)(n - 1) 



(C.7) 



(a, (3 - 1 - A, a - 1, j3 - A, a) 



(a - - A,a,/3 - A, a) 



r(/* + i) 



7r^a(a - l)a(a)a(/3) 



a — 1 



(C.8) 



r(M) 



- + £(/?) -3 



- B(a-1) + 



/? - 1 _ 1" 
a — 1 /Lt. 



(C.9) 



(a,/3 - A,a,/3 - 1 - A, a) 



TT 2 ^a 2 (a)a((3 - 1) 
1 «(/?-!) 



+ £(/?) - B(a - 1) 



2(a-l)(ji-0) 
(3/i-l)(a-l) 2" 
2fi(fi-P) n_ 



(CIO) 



(a - 1, /? + 1 - A, a, /? - A, a - 1} 



vr 2 ^a 2 (a-l)a(/3 + l) 
(/x - a)^r(// - 1) 
+ 1 /?( / u-/3 + a~2) - 



(a-1) 



(C.ll) 
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{a - 1, (3 - 1 - A, a, (3 - A, a) 



Tx 2 ^a{a - l)a(a)a((3) 
ATM 

+ A (b(0) - B(a - 1) 



09- 1) ' 
(a-1). 



(C.12) 



(a,/3 - 2 - A, a - - A, a) = 



tt 2 ^(2 - n)a(a - l)a(a)a(f3) 
IV+T) 



+ ^ 2 + /U-l , 2(/i-/3) 
(a-1) 



(C.13) 



(a,/3 - 2 - A,a,/3 - A, a) 



+ 



2n 2 ^a(a - l)a(a)a(/3) 
ir 2 »a(a)a(a - l)a(p) 
1 4(/i - 0) 



(/3- 1)^-1) 



M(« - 1) 2 



(a - 1) (a - 1) 

^ 2 + I s 



- B(a-1) + 



n(n - 1) 



(C.14) 



(a + 1, /3 - 2 - A, a, f3 - 1 - A, a) 



x 



7r 2 ^a 2 (a)a(/3-l)(/3-l) 
2(/z - a - l)T(fi + 1) 
(0-2)09-3) 



2a 

4(a- !)(/?- 2) 
09-1) 



+ /X-/3+1 

(C.15) 



(a + l,/3 - 1 - A, a, (3 - A, a - 1} = 



a)a(a — l)a(/3 — 1) 
(/x - a - l)r(/i + 1) 

"/U — Q + 1 (/i — 1) 



(3-1 



a 



(C.16) 
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(a + 1, (3 - 1 - A, a, /3 - 1 - A, a) 
_ n 2 ^a(a)a(a - l)a((3 - 1) 
r(//)(a- 1)03-1) 

+ 2(B(J3)-B(a-1)) + 



4 2 
A + /3-1 



2 2 
-- + - 

a fj, — a — 1 



(/?-!) 



-/?)(«-!) 
4(2a + 2 j 9-3)(a-/A-l)(a-l) | 



l){n-a- 1) 



4(/x + a - 3) 
(/i - 1)0 - a - 1) 



2(a + /9 - l)(a - // - !)(/? - 2) f/?-3 + /z-/?+l 



//(//- l)(n-a- 1) 
(/3-l) 2 0-/3)0-a + l) 



— 1)0 — a — 1) V/9 — 1 (/i — a + 1) 



a 



2a /3 - 2 



(C.17) 
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Figure Captions. 

Fig. 1. Skeleton Dyson equation for the electron. 

Fig. 2. Skeleton Dyson equation for the photon. 

Fig. 3. Uniqueness rule for a bosonic vertex. 

Fig. 4. Electron photon vertex. 

Fig. 5. Integration by parts rule for gauge vertex. 

Fig. 6. Additional integration by parts rule. 

Fig. 7. Further integration by parts rule for gauge vertex. 

Fig. 8. Graphical representation of IL^ (x) . 

Fig. 9. Photon self energy after integrating by parts. 

Fig. 10. Basic two loop self energy graph. 

Fig. 11. Various basic bosonic graphs contributing to 

Fig. 12. Basic rule for recursion relations. 
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